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The quantum theory of the fluctuations of the van der Waals (vdW) force between 
macroscopic bodies is developed. Unlike the mean vdW force that is determined by all 
quantum states that contribute to the optical absorption, the energies of those excitations 
of the interacting bodies that contribute effectively to the vdW force noise are thermal 
energies. Contrary to the mean vdW force the vdW force noise drops with decreasing 
temperature. Due to these differences the main mechanism of the mean vdW force and 
that of the vdW force noise may be different, e.g., the mean vdW force is determined 
by electronic excitations, the force noise by the random lattice or impurity dynamics. 
Since the vdW force is linear in the fields squared the dependence of the vdW force noise 
correlation function on time difference is determined not by one but by two frequencies: the 
spectral density depends not only on the frequency of its measurement us but also on the 
integral over the second frequency. The dependence of its integrand on this frequency and 
temperature is quantum-like even at small us <C ksT/h, and such may be the behaviour 
of the vdW noise. In the case of interacting metals the vdW force noise is enhanced by dc 
voltage applied across these metals. This additional noise is proportional to the voltage 
squared and to the spectral density of the random electric field at the frequency of noise 
measurement. The theory is in qualitative agreement with experiments. 

PACS numbers: 05.40.Ca, 07.79.Lh, 72.70.+m 

1. INTRODUCTION 

After the fundamental works on Brownian motion by A. Einstein and M. Smolu- 
chowski in the early 20-th century and after the experiments and theory on electrical noise 
in 1920s by J.B. Johnson and H. Nyquist it was realized that the ultimate limit of the 
sensitivity of any device and the accuracy of any measurement is determined by the fluc- 
tuations (noise) in the main physical phenomenon on which the device or measurement is 
based (not, for instance, by noise in the auxiliary circuitry). In the last two decades sev- 
eral types of scanning microscopy with atomic resolution have been developed: Scanning 
Tunneling Microscopy (STM), Scanning Atomic Force Microscopy (SAFM), and Magnetic 
Resonance Force Microscopy (MRFM). The Scanning Force Microscopies are based on the 
measurement of the force between the sample and the tip (probe) mounted on the end of a 
cantilever using a sensitive laser interferometry that measures the bending of the cantilever 
under the force. At present, owing to the progress of experimental methods, it is possible 
to measure very small forces on the order or even less than 10 -18 N. 1-3 Recently the force 
noise in devices typical for SAFM have been measured. ' 5 The spectral density of the force 
noise was found to be ~ 10 -35 — 10 -32 N 2 /Hz. The spectral density of the noise depends 
on the sample material (Ref. 5) and increases as the vacuum gap L between the probe 
and the sample is decreased. 4 ' 5 There are several questions that arise from these studies. 



1 



1. What is the minimum force measurable by nanomicroscopies? How does this 
ultimate limit depend on the materials, dimensions, and geometry of the tip and sample? 

2. What is the temperature dependence of the van der Waals force noise? 

3. Is the mechanism of the measured force noise inherent in the very mechanism of the 
measured force (and cannot be eliminated) or is it caused by extraneous mechanisms that 
can be eliminated by cleaning the surface of the sample, improving the cantilever quality 
and so on? 

The main force between neutral bodies at distances between them in the range ~ 1 
nm — 1 fxm is the van der Waals (below, for brevity, vdW) force. The vdW force between 
neutral atoms and molecules 6 , atoms and lossless metals 7 or macroscopic bodies 8 is the 
result of correlation between the motion of charged particles in the interacting bodies. 
This correlation originates due to interaction via electromagnetic fields. In macroscopic 
bodies these fields are generated by random current (polarization) sources. E.M. Lifshitz 
calculated the mean vdW force between macroscopic bodies (per unit area) as the mean 
value of the Maxwell stress tensor, T^-, in the vacuum gap near the surface of one of the 
interacting bodies (Ref. 8). The sources of Maxwell surface stress are all charges and 
currents inside the body enveloped by this surface: 

T ij(^t) = ±[(-E i E j + l -5 l3 E 2 ) + (-HiHj + \^H 2 )\ (1.1) 

Here and below r is the 2D radius vector parallel to the surface, E(r, t) and H(r, t) are the 
randomly varying electric and magnetic fields at the surface, respectively. In the first 
index i corresponds to the direction of the force, the second to the direction perpendicular 
to the element of the surface. In a realistic device this Maxwell stress should be calculated 
for the surface of the tip (in Ref. 8 a mathematically simpler model was used). At a given 
distance L between the probe's point and the sample surface (perpendicular to the OX 
axis), the force acting on the probe in the OX direction is an integral of the Maxwell stress 
T x j over the tip's surface, A: 

F x (t) = [ T XJ (t)dA J: (1.2) 



J A 

where dA is an element of the surface. For the system of two parallel semiinfinite slabs 
separated by a vacuum gap L, at temperature T — > and small L, the mean vdW force 
per unit surface area reads 8 : 
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Here £ is a real frequency variable. The dielectric permittivity at imaginary frequency, 
e(i£), is a real quantity that monotonously falls off with £ from the value of the static 
dielectric permittivity e(0) down to e(zoo) = 1, integration over dimensionless p originates 
from the integration over the wave vector of the fields. Equation (1.3) is valid at small 



2 



distances L between the bodies when the retardation time of the electromagnetic waves 
traveling between the bodies can be neglected. 8 It happens when L <C c/ui where c is the 
light speed and u> is the maximum frequency of significant optical absorption. Practically, 
hu ~ 10 - 100 eV. 

As follows from the theory 6 ' 8 , the mean vdW force between atoms or macroscopic 
bodies is determined by the entire spectrum of their optical absorption. The main part of 
this wide spectrum usually corresponds to frequencies u ^> ksT/h. At these frequencies, 
the quantum states that are the final states of the optical transitions are not populated. 
Therefore the variation of the temperature does not affect significantly the mean vdW 
force. We show below that the range of excitations' energies that determines the vdW 
force noise and its temperature dependence is drastically different from the energy range 
that determines the mean vdW force. 

The force noise is measured by observing the random displacements of the cantilever. 
There are two mechanisms of random cantilever motion: 1) random motion of the me- 
chanical system due to internal friction (damping) affected by the vdW interaction with 
the sample and 2) random variation of the electromagnetic fields that are the sources of 
the vdW force. The results obtained in the first problem are presented and analyzed in 
the comprehensive review by Giessibl 9 (see also 3 ). This article is devoted to the second 
problem, i.e., to the theory of "purely" vdW noise. 

2. QUANTUM CORRELATION FUNCTIONS 

Let F(t) be a fluctuating physical quantity, e.g. force, and (F) its mean value. Then 
SF(t) = F(t) — (F) is the time-dependent fluctuation of the quantity and C(ti,t 2 ) = 
(SF(ti) 5F(t2)) is the correlation function of the fluctuations (the angle brackets denote the 
statistical and quantum averaging). In general, F(t) is a Heisenberg quantum-mechanical 
operator, F{t). Two such operators, 8F{t\) and 8F(t 2 ), taken at different times t\ ^ t 2 (as 
in the correlation function), do not commute and C(ti,t 2 ) depends on their order. There 
are actually two correlation functions 10-12 : 

C-(t u t 2 ) = C_(ti - t 2 ) = {SFitJSFfo)) = (SFfa -t 2 + t )8F(t )), 

C+(t lt t 2 ) ee C+(h - t 2 ) = (5F(t 2 )5F(t 1 )) = (5F(t )5F(t 1 -t 2 + t )) = C_ {-(h - t 2 )}. 

(2-1) 

According to the Wiener-Khintchine theorem, the spectral density is twice the Fourier 
transform of the correlation function. It follows that actually there are two Fourier trans- 
forms of the correlation functions and two spectral densities: 

/+00 
dte iult C±(t), S±(f) = 2C±(oj) t = tx-t 2 . (2.2) 
-00 

Here / = uj/2tx is the frequency assumed positive. Since C_(£) = C + (— t), one can also 
present S-(f) and C-(u) as the functions S + (f) and C + (uj) at the frequency with opposite 
sign: 

5_(/) = 5+(-/), C_{u) = C + {-u). (2.2a) 
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The correlation functions and spectral densities can be presented in terms of the 
eigenstates \m) of the system, their energies E m , and the states' probabilities w m which 
for equilibrium systems are given by the Gibbs distribution: 

C + (t) = Y t w m (rn\SF(0)6F(t)\m) 

m 

C_(£) = J2 w m(m\5F(t)5F(0)\m) 

m 

S±(f) = 2hJ2^m\(SF) mn \ 2 5(E rn 

mn 

In an equilibrium system w m oc exp[— E m /kBT] and 

S-(f)=^P[^/k B T]S+(f). (2.4) 

The spectral density that corresponds to the commonly used symmetrized correlation 
function equals S s (f) = ^[S + (f) + S-(f)]. Using Eq. (2.4) one can easily express S± in 
terms of S s and vice versa. 

The arguments of the delta functions in Eq. (2.3) show that S+(f) corresponds to 
the emission of energy from the sample (system) under study and S- (/) to the absorption 
of energy. 10-12 Each spectral density is proportional to the rate of emission or absorption 
under the condition that there is no backward reflection of the energy. The actually 
measured spectral density depends on the properties of the measuring device and the 
measurement conditions. If the sample under study is at absolute zero temperature it 
cannot radiate energy or momentum to the measuring device and the spectral density 
S + (f) = 0. However, if at the same time the device temperature > the sample 
absorbs the noise energy emitted by the device, S-(f) ^ 0. In principle, this spectrum 
also may be used as a characteristic noise spectrum of the sample. 10-12 If Td = the 
measured "absorption" noise is also zero. 

The spectral density S + (f) of the vdW force noise is zero at T = unlike the mean 
vdW force (F) that is non-zero at the same zero temperature (see Eq. (1.3)). This 
conclusion is in full agreement with the fact that at T = the system, being in the ground 
state, cannot transfer energy or momentum to the noise detector (see above). It means 
that the vdW force noise grows with temperature T. If the experiment actually measures 
S+(f) its value at T = is zero. The quantum theory of the force noise is instructive for 
understanding the fundamental difference between the properties of the mean vdW force 
and those of the vdW force noise. 

3. THE VAN DER WAALS FORCE NOISE 

IN THE SYSTEM OF TWO SEMIINFINITE SLABS 

^From the mathematical point of view, the simplest model of a system of macroscopic 
bodies interacting via vdW force is the system of two semiinfinite slabs with parallel 
surfaces separated by a vacuum gap. It was introduced by E.M. Lifshitz 8 and later used 
by many authors. We shall use this model below to calculate the spectral densities S±(f) 



= J2 w rn\(6F) mn fe^- E ^, 

mn 

= J2 w m\(SF) mn \^ E -- E ^\ (2 . 3) 

mn 

-E nT hf). 
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of fluctuating vdW surface stress and force. It is convenient in this model to introduce 
the 2D radius-vector r in the surface plane of the slabs and the OX axis normal to this 
surface. The Langevin source in this problem is K(x, r, t) = 47rP(x, r, t), where P is the 
extraneous (source) polarization vector. In Ref. 8 the Fourier expansion of P in both slabs 
is introduced: 

AnP m {x, r,t)= [ rfqe iqr f°° dk x { cos (k x x)g m (k q t), if n = 1; 

J J x \co8[k x (x-L)]g m (k x ,<i,t), ifn = 2. v 1 

g m (k x ,ci,t) = J dre" !qr J dxX n (k x ; x)P m (x, r, t), (3.2) 

where 

dxX n {k x ;x)...= i I < Lo dxcos ( k - x )--- ifn = l; 

\ dx cos[k x (x — L)] . . . if n = 2. 

In Eqs. (3.1) and (3.2) q is the wave vector in the plane of the surface, % and n number the 
vector components and the slabs, respectively, m = (i,n), g m (k x , q, t) is the cosine Fourier 
transform of the Langevin source in the n-th slab. As follows from Eq. (3.2) the Fourier 
amplitudes g n (k x , q, u) = g n (—k x , q, uj), i-e., they are even functions of k x . 

We are using below the version of the Langevin approach used by Lifshitz 8 . The elec- 
tric field, E(x, r, t), magnetic field H(x, r, t), the stresses Ty(r,£) at the surfaces, and the 
Langevin sources of the fields are considered as quantum-mechanical operators satisfying 
the same equations as their classical counterparts (a comprehensive review of this approach 
is given by Henry and Kazarinov 13 ). 

The correlation functions of the g n sources can be expressed in terms of the po- 
larization correlators. In general, the Fourier transform of the polarization correlator is 
proportional to the imaginary part of the dielectric permittivity €ij(xi, ri; x 2 , r 2 ; u;) in the 
linear integral relationship between the electric displacement Di(xi, r±, uj) and the electric 
field Ej{x2,Y2,uj) (see 14 ). In the macroscopic approximation, 

tij{xi, r x ; x 2 , r 2 ; uj) 8{x\ - x 2 )8{r 1 - r 2 )e ij -(w). (3.4) 

Due to 8{x\ — X2) the double integral J dxiX n (k x i;xi) J cte 2 X n (fc x2 ; X2) in the correlator 
(gm 1 {k x i,c[ 1 ,uj 1 )g m2 (k X 2,q.2,uj2)) (Eq. (3.2))is reduced to one integral: 



/ 



7T 

dxX n (k xl ;x)X n (k x2 ; x) = -8(k xl - k x2 ). (3.5) 



As follows from Eqs. (2. 2), (3. 2), (3.4), and (3.5) the correlators of the random quan- 
tities g m (k x ,q,uj) are: 

(9m 1 (k x i,qi,uj 1 )g m2 (k x2 ,q2,uJ2)) = 2 4 5(k xl - /c x2 )5(^i + cu 2 )5(qi + q 2 )C_ mi (u>i) 
(9m 1 ( k xi,<li,u 1 )g? n2 (k X 2 1 q 2 ,uj2)) = 2 4 5(k xl - k x i)8{u\ + w 2 )5(qi + q 2 )C_ mi (-u>i), 
(9m 1 ( k xi,<li,uj 1 )g rn2 (k X 2,q.2,UJ2)) = 2 4 5(k xl - k X 2)6(ui - £J 2 )5(qi - q 2 )C_ mi (-cui), 
(g mi {k x i,q. 1 ,uj 1 )g^ n2 (k X 2,q.2,UJ2)) = 2 4 5(k xl - k X 2)6(ui - £J 2 )5(qi - q 2 )C_ mi (cui). 

(3.6) 
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For brevity we omitted 5 mim2 in the rhs of these equations. Equations (3.6) are presented 
for positive frequency arguments of the functions C_ . If the frequency argument is negative 
this function is in fact a C+ function (Eq. (2.2a)). The equilbrium correlation function of 
polarization is well known (see, for instance, 14 and Eq. (3.4)). For an isotropic solid, 

C T , m (u) = J d( Xl - x 2 )d(r 1 - v 2 )d{h - t 2 )e ±iuJ ^-^ (P m (x u r 1} h)P m {x 2 , r 2 , t 2 )) 

_t^„ (N(uj) + 1 forC_, m (u;); 
" 2n en[UJ) \N(uj) for C+, m (oj). ' 

(3-7) 

Here N(u) is the Planck equation for the number of photons in one mode. 

The vdW force is determined by the component T xx (x = 0, r, t) of the Maxwell stress 
tensor: 

T xx (r, t) = ± [El + E 2 Z -E* + H 2 y + H 2 Z - H 2 X ] , (3.8) 

where the fields are taken at x = 0. The Fourier components of these fields are calculated 
in Ref. 8. At each given q: 

E x (q, to) = -(q/p)(v y - w y ), E y (q, u) = v y + w y , E z (q, u) = v z + w z , 

c c to 1 

H x (q,(j) = -q(y z + w z ), H y (q,uj) = p(v z -w z ), H z (q,w) = --(v y - w y ). 

(3-9) 

Here p = a/cu 2 /c 2 — q 2 , the vectors v and w depend on the wave vector q and the frequency 
to. They are integrals of the Langevin sources. 8 Introducing the notation 

/+oo 
dk x g m (k x ,q,w)/(kl- s*), (3.10) 
-oo 

where s 2 = (uj/c) 2 e n (uj) — q 2 , we present them as follows: 

v y + w y = \e~ wL (e 2 p + s 2 ) - e ipL {e 2 p - s 2 ) [qg lx - s\g\ v ] - 2s!S 2 [qg 2x + s 2 g 2y ]^, 

v y -w y = ^-{si \e~ wL (e 2 p + s 2 ) + e ipL (e 2 p - s 2 ) [qg lx - stg ly \ + 2s 2 e 1 p[qg 2x + s 2 g 2y ]}. 



,2 



v z + w z = -2p-^— [sie ipL g lz + s 2 g 2z ], v z - w z = 2-^— sis 2 [-e ipL gi z + hz], 
c A c A 

A = e ipL ( Sl - eip)(s 2 - e 2 p) - e" ipL (si + e lP )(s 2 + e 2 p), 
A' = e * pL ( Sl -p)(s 2 -p) - e- ipL ( Sl +p)(s 2 +p). 

(3.H) 

The vectors v and w are linear in g m (q : u). The coefficients do not depend on k x . 
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As follows from Eqs. (3.8) and (3.9) the unaveraged stress at the surface x = equals: 
T xx (r,t) = ±-J |^*i*iVK^^ 

+ [(ujuj'/c 2 - qq')/pp'*}(v y - w y )(q,u)(v y - w y )* (q_' , u') 
+ [l-(c 2 qq , /ojuj , )}(v z + w z )(v z + w z )*(c i , ,oj') 

+ (c 2 pp'*/uu')(v z -w z )(q,u)(v z -w z y(q',u')y 

(3.12) 

After substitution of Eqs. (3.11) in Eq. (3.12) we obtain the unaveraged stress 
T xx (r, t) as a sum of terms each of which is proportional to a product # m (q, cu)^j' re ;(q / , ui'). 
The equation for T xx takes the form: 

T xx (r,t) = i- / ^^ eiKQ-QOr-^-.O*] £ aw (q, q', c')^(q, -)^(q', c')- 

m,m' 

(3.13) 

The coefficients a mm / are calculated in the Appendix using Eqs. (3.11) and (3.12). Since 
the stress must be real (Hermitian) these coefficients must satisfy the relationship: 

a mm /(q,w;q',c/) = a* m , m {q[, J\ q, to). (3.14) 

The unaveraged T xx (Eq. (3.12)) depends on uj and q in two ways: either directly on 
u 2 and q 2 or on arguments of the dielectric permittivities, e n (q, a;). It depends on the 
signs of q and u> only in the second way. Since e n (— q, —uj) = e* (q, a;) we obtain one more 
relationship: 

a mm /(-q, -u;;-q',-u/) = c4 m ,(q,w;q',c/). (3.15) 

Our goal is to calculate the correlation functions C±^T xx T xx (^i,ti;r2,t2) and the cor- 
responding spectral densities (Eq. (2.2)) of the surface stress fluctuations. Due to general 
relationship between the spectral densities S-(f) and S+(f) (Eq. (2.4)) the problem may 
be reduced to the calculation of C + and the corresponding S + (f): 

C +! T xxTxx (r 1 ,t 1 ;r 2 ,t 2 ) = (T x ,(r 2 , ^T^, tx)> - (T xx (0, 0)> (T xx (0, 0)>. (3.16) 

Each correlation function is expressed in terms of the mean product of four random 
Langevin sources g m . The Langevin sources in a macroscopic body are usually considered 
Gaussian. It follows that the mean value of this product can be decoupled and presented 
as a sum of products of pair correlation functions with the order of operators g m preserved 
within each pair. The correlation function equals: 

dqi dq^ dc\2 dq 2 <ko\(k>j' x du>2cLu 2 



C+,T xx T xx {ruh\ r 2 ,£ 2 ) - {^) J 



e 



(2tt)< 

:[(qi-q' 1 )ri + (q 2 -q2)r2-(wi-Wi) t l-(^2-W2)t 2 ] 



X 5^ am i m i( qi ' Ci;i;q/ l' Ci; l) am 2mi ! (q2,W2;q2,W2) ( 3 - 17 ) 
X [(<7m 2 (q2, W 2 )^(q 2 , ^ 2 )#mi ( k xl, <ll, ^iWm'^ fai » 
- {9m 2 {(\2^2)g*m'S^2^2)) (drm (qi , ^l)9 m ' (Qi » "l)) 



After decoupling, the expression in the rectangular brackets takes the form: 

(q2, w 2 )5 mi (qi, wi)> {g* m ' (q^km' (qWi)) 

/ / ( 3 - 18 ) 

+ 5 m ' im2 5 mim ^ (^ 1 (q 2 ,^2)^m 1 (qi,^i)) (5m' 1 (q2»W2)^' i (qi,o; 1 )). 

Using Eq. (3.6) the correlators of the g n functions are expressed in terms of the 
correlators of the polarization sources of fluctuations: 



(<7m(q2, <*>2)<7m(qi, = 24< ^(qi + Q2)5(^l + k>2)C+,m(k>l)-7n(?l, ^l), 

(5m(q2» w 2)5m(qi» w i)) = 2 4 Wi + Oa^C^i + u; 2 )C_ m (u>i)J n (^,^i), 
(^m(q2. ^2)^m(qi,^i)) = 2 4 5(qi - q 2 )5(^i - W2)C + , m (a;i)/ n (gi, Wi), 

(5m (02, w 2 )^(qi, wj)) = 2 4 <Kq' 1 - q2)^K - w 2 )C- mKKn^i, wj). 



(3.19) 



Here 



If - /" + °° dkx - 71 ^ (3 20) 

J-oo \ k l - sl(q,u)][kl - sl(q, -oj)} \s n (q, uj)\ 2 Ims n (q, u) ' 

Thus the surface stress correlation function equals: 
x a mm /(q,w;q',w') [a mm /(-q, -a;; -q', -a/) + a m / m (q', u/; q, w)] (3.21) 

mm' 

x C+,m(q, w)C_, m '(-q', Lv')I n (q, u)I n i(q', u') 
According to Eqs. (3.14) and (3.15) the correlation function takes the form: 

C + ,r«r„(r lltl ; r2 , i2 ) = A/ ^^JI^e'Ko-'X-- )-<— ' ><«-«l 
x |a mm /(q,w;q ,u')\ 2 C + , m {u)C-, m >{u')I n {q,u)I n >{q\u'). 



(3.22) 



The spectral density is twice the Fourier transform of the correlation function. Let 
us be the frequency at which the spectral density is measured. After the substitution 
u> = O + W5/2, ui' = O — ujs/2 the integration over frequencies is reduced to the integration 
over O: 

S +,t xx t xx (ri - r 2 ,u) S ) 

= ^4 / 7^ / ^qV^')( ri - r2 ) Yl |amm'(q,0 + ^/2;q',0-^/2)| 2 ( 3 .23) 

— 00 mm' 

x C+, m (0 + cj s /2)C_ im ,(fi - u s /2)I n {q, O + to s /2)I nt (q', O - u; 5 /2). 



The only temperature dependent quantities in Eq. (3.23) are the correlation functions C±. 

Using the relation C±, m (-u;) = C Ttm (u) (Eq. (2.2a)) and Eqs. (3.14), (3.15), (3.20) 
it is easy to show that the integrand remains invariant under substitutions O — > — O, 
q < — > q'. Therefore we can integrate over positive frequencies O only: 

x \amm>(ci,n + uJs/2;ci , ,n-us/2)\ 2 I n (q, Q + u s /2)I n ,(q' ,Q - u s /2) (3-24) 

mm' 

x C+, m (0 + u s /2)C-, m/ (n - us/2). 

In experiments, not the stress but the force acting onto an area of the body is mea- 
sured: F(t) = f drT xx (r,t) (see Eq. (1.2)). T xx depends on r through the exponential 
function e^ q_q ) r in the integrand. For a circular area of radius r c the integration over r 
replaces this function with b = 27rr c Ji(|q — q'|r c )/|q — q'|. The spectral density of this 
force, Sf(ws), is given by Eq. (3.24) in which the exponential function of ri — r2 is replaced 
by b 2 : 

Sf(oj S ) = ^J o ^ J dqdq'Mxdq- q'|r c )/|q- q'|] 2 

x l a — '(q, n + ws/2;q'fi - uj s /2)\ 2 I n (q, Q + u s /2)I n ,(q' , Sl-u s /2) (3-25) 

mm' 

x C+, m (0 + u s /2)C-, m ,(n - u s /2). 

The integration in Eqs. (3.24) and (3.25) is simplified owing to the dependence of the 
correlation functions only on O while other factors in the integrand, at O/c L _1 , depend, 
in fact, only on q and q' . 

The integrand in the equations for the spectral densities of T xx and vdW force is 
proportional to a product of two Langevin sources' spectral densities: one of C + type with 
frequency argument O + us/2, the other of C_ type with frequency O — us/2. In that 
part of the integrand where O < us/2 the frequency argument of C_ is negative, i.e., 
according to Eq. (2.2a), both functions are of C + type and both tend to zero as T — > 
(Eq. (2.3)). However in known experiments the frequency us at which the spectral density 
is measured (~ 10 4 — 10 5 Hz) is much smaller than the frequencies of strong dissipation 
e'n(ui) because, in general, e"(0) = 0. The main contribution to the force noise comes from 
O ^> us- Then the temperature dependence of the vdW spectral density is determined 
by the product C + (0)C_(0) oc N(Q)[N(Q) + 1]. Due to N(Q), the contribution of the 
dissipation e"(0) at O > k B T/% is exponentially small and does not contribute to the 
vdW force noise despite it fully contributes to the mean vdW force. Thus each term in 
the integrand of S^t^t^ falls off to zero as the temperature T — > unlike the mean 
vdW force that tends to a non-zero value (Eq. (1.3)). An interesting peculiarity of the 
vdW force noise is also worth mentioning: even at small frequencies us ksT/h, i.e. at 
conditions that are commonly considered as classical, the temperature dependence of the 
noise, due to the function C + (0 + a;s/2) in the integrand, behaves, at O > kBT/%, as that 
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of a quantum noise. However, if at us < <~ k B T/h the dissipation e^(O) is significant 
the temperature dependence of the force noise is ~ T 2 . 

The analysis presented above shows that only the excitations with energies Tiui smaller 
or on the order of thermal energy k B T effectively contribute to the vdW force noise. In 
this respect the vdW force noise differs from the mean vdW force which is determined 
by all excitations that contribute significantly to the electromagnetic dissipation in the 
interacting bodies, regardless of their energy in units k B T (Eq. (1.3)). Due to the fact 
that the contribution of the thermal excitation energies ~ k B T to the mean vdW force is 
usually insignificant and, on the contrary, only these energies contribute significantly to 
the vdW force noise the mechanisms of these two vdW phenomena may be quite different. 
The mechanism of the mean vdW force is mainly connected with electronic transitions, 
that of vdW force noise may result mainly from lattice vibrations, e.g., transverse optical 
phonons, from the motion of impurities or from quasi-stationary currents in metals. The 
two measured quantities, the mean force and the force noise, supplement each other. 

Other factors in Eqs. (3.24) and (3.25) (besides C± iTO ) depend on fi/c. When fi/c <C 
L -1 the ratio fi/c can be dropped and these factors depend on q and q' only. Then (see 
Ref. 8), 

p » iq; s n » p; e ±ipL = e^ L ; I n {q, w)I n ,trf , u') oc {qq')-\ (3-26) 

The elements of the matrix a mm '(q, oj; q',a/) are calculated in the Appendix. In the 
approximation of small (OL/c) 2 , the elements |a n i )TJ ,v| 2 with = x,y are proportional 
to (qq') 4 and a ZUjZn ' oc (O/c) 2 can be dropped. Then the part of the integrand in Eq. 
(3.25) that depends on q and q' can be presented as a function of qL and q'L. One can 
see that if r c ~ L the spectral density of the vdW force is oc L~ 2 . 

The magnitude of the vdW force noise is determined, like that of the mean vdW 
force, by the complex dielectric permittivities of the interacting bodies. As is well known, 
in heteropolar dielectric crystals a strong IR absorption peak at the frequency of the 
transverse optical phonons, u t , exists. In the range of this peak: 

(f'Mn^utieo-eooMu-ut)- (3-27) 

The frequency ui t is ~ 10 14 s _1 , i.e., Tiu t ~ ksT at room temperature T. This dissipation is 
a potential mechanism of significant vdW force noise. Substituting [e"(0)] 2 in Eq. (3.25) 
we must replace the 5-function peak with some more realistic, e.g. Lorentzian peak of 
finite width 7. In Eq. (3.25) and preceeding equations, the integral over the peak squared 
yields: 

7RJ 2 (e - Coo) 2 



/ 



dQ[e"(n)Y * , (3.28) 

Then the estimate of the force spectral density at small lvs is: 

q , x [Hut (eo ~ Coo)] 2 ) 1 1 ( ooq\ 

S) 7L 2 e W [(^t/k B T) -11- eM(^t/k B TY { ' 1 

Substitution of characteristic values, eo — ~ 1, L ~ 1 nm, 7 ~ 0.1a; t , yields: 

1 IN 2 

S±AUS) ~ eM(^ t /k B T) -11- e X p[-(^ t /k B T) ' (3 - 30) 
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4. EFFECT OF dc VOLTAGE BETWEEN MACROSCOPIC METALS 
ON THE VAN DER WAALS FORCE NOISE 

Consider a system of a probe and a sample that are both metallic or metallized. In 
this case a dc voltage U can be applied across the vacuum gap between the probe and 
the sample. The goal is to find the effect of the dc voltage on the mean force and on the 
force noise. The mean stress (and force) is simply a sum of vdW stress (Tij) vdW and the 
electrostatic stress due to the applied dc voltage, (Ty)u- The latter is proportional to U 2 . 
These two forces are simply additive. However, the vdW force noise is strongly affected 
by the application of the dc voltage. 5 

Let us calculate the change of the stress correlation function with voltage. The total 
electric field at any point r in vacuum is a sum of the dc field F(r) and the randomly 
fluctuating field E(r, t): 

E tot (r,t)=F(r)+E(r,t). (4.1) 
The total stress is a sum of three terms: 

Tij (r, t) = T$> (v,t) + T^\r,t) + (r, t) . (4.2) 

Here the upper index is the total power of random field components in the term. As follows 
from Eq. (1.1), 

(2) 

and Tj ■ is the stress due to the random electric and magnetic fields only. The correlation 
function of the total stress, 

<7 Ti ,T fcI (ri,*i;r 2 ,* 2 ) = (^(n, £i)T fc/ (r 2 , t 2 )) - Cfy(n,*i)) (2fc/(r 2 ,£ 2 )). (4.4) 

We substitute Eq. (4.2) in this equation and take into account that in equilibrium the 
mean value of products of odd number of random fields is zero. The variation of the 
correlation function due to the static electric field (dc voltage U) equals: 

4^ T >i,ti;r 2 ,t 2 ) = (T^^M)^ (T2M)) 

= (4^((^' F(ri) ' E(ri ^ l} " F *( r O s i( r i»*0 " fi(ri)£ 4 (n,fi)) (4-5) 
x (<5 fcZ F(r 2 ) -E(r 2 ,£ 2 ) - F fc (r 2 )^(r 2 , £ 2 ) - F z (r 2 )E fc (r 2 , t 2 ))}. 

Obviously, at any geometry of the probe and sample, the spectral density of the additional 
noise due to dc voltage, S^ u \f) : is proportional to U 2 , does not depend on the sign of 
the voltage, and is linear in the spectral densities of the random electric field E compo- 
nents, SEiEj (/; r i — r 2), at the frequency of the spectral density of noise measurement (no 
additional integration over frequency). It means that the electric field spectral densities 
determine the steepness of the additional noise growth with U 2 . At "classical" frequencies 
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/ <C k B T/h the spectral densities SsiE-if) are roughly linear in temperature T. Under 
these conditions the noise (f) grows linearly with temperature. It must be emphasized 
that the noise oc U 2 is created by a noiseless dc voltage and therefore can be viewed as an 
amplification of noise generated by the only its sources, i.e., by the sources of the random 
electromagnetic fields. 

In the model of two semiinfinite parallel metals separated by a vacuum gap of thickness 
L the dc field in vacuum F x = (U/L), and the vdW force is determined only by T xx . It 
follows from Eq. (4.5) that the spectral density of the stress equals: 

Ck,(/;ri - r 2 ) = ^S EmE .{f;ri ~ (4.6) 

The problem has been reduced to the calculation of the spectral density of the electric field 
E x in the vacuum gap near the surface of any of the two slabs. 

5. MAIN CONCLUSIONS AND COMPARISON WITH EXPERIMENTS 

Since the geometry of the theoretical model (two semiinfinite slabs) strongly differs 
from that of the real SAFM device (acute metallic tip above the sample) the quantitative 
comparison of the theory and experiment is not justified. Therefore we present below a 
summary of qualitative results and compare them with the known experimental data. 

First, we found that the properties of the mean van der Waals force and the wdW 
force noise strongly differ. The former are determined by all states of the interacting bod- 
ies that contribute to the optical transitions in these bodies. At common temperatures of 
measurement, the major part of these quantum states are not populated. The mean vdW 
force is therefore non-zero at zero absolute temperature and the variation of the tempera- 
ture only weakly affects this zero-temperature value. On the contrary, the measurement of 
the force noise requires transfer of energy and / or momentum from the noisy system to the 
noise detector. Only the quantum states with energies (reckoned from the ground state) 
less or on the order of k B T contribute to the force noise. That is why the force noise, 
unlike the mean force, falls off at T — > 0. At low temperatures it falls off exponentially. 
This conclusion is qualitatively supported by experimental findings 5 : in the system Au 
film— Au-coated cantilever the effective damping coefficient V = Sp/^k B T was found to 
be ~6 times lower at 77 K than at 300 K. It means that the spectral density of the force 
noise was ~24 times smaller at temperature only 4 times lower. 

Van der Waals forces between two bodies are due to the correlation of charges and 
currents in both bodies via electromagnetic fields. Since the system's excitations with 
energies <~ k B T comprise a very small contribution to the mean vdW force and, on the 
contrary, they comprise the main contribution to the vdW force noise, the mechanisms of 
the electromagnetic fields' sources, i.e. the origin of the vdW forces may be quite different 
in both cases. 

The vdW force between macroscopic bodies is expressed in terms of the Maxwell 
surface stress, Tij(r, t), that is bilinear in the components of the randomly varying electric 
and magnetic fields. Therefore the spectral density of the vdW force noise depends on 
the (given) frequency us at which it is measured but is also an integral over the second 
frequency, O (Sec. 3). The integrand is proportional to the product of two spectral densities 
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of the random electromagnetic fields' Langevin sources. In the range of small wg, much 
smaller than the frequencies of significant dissipation (e" 2 ), the spectral density of the 
vdW force noise is practically independent of us- Then the temperature dependence of 
the force noise is determined by the integral over Q the integrand of which depends on 
exp[hQ/kBT]. If the frequencies of strong dissipation are > ksT/Ji the noise is quantum 
even at us <C ksT/h. 

In experiments the dependence of the vdW noise on L is ~ L~ d with d = 3 (Ref. 4) 
and d = 1.3 (Ref. 5), respectively (the difference is due to different geometry of the used 
systems "cantilever— sample" 5 ). According to Sec. 3, in the two-slab model the force noise 
falls off with L also as ~ L~ d but with d ~ 2. This dependence is close to the experimental 
ones. 

The striking phenomenon of force noise enhancement between two metals by dc 
voltage 5 follows immediately from the general equation for the vdW force (Sec. 4). Ac- 
cording to the experiment 5 and to the theory the additional force noise is proportional to 
the dc voltage squared. This agreement between the theory and experiment proves once 
more that both the force and force noise between the two metals are of vdW origin. 

As is shown in Sec. 3, the strong IR absorption peak at the frequency of the trans- 
verse optical phonons, uit, may be a significant source of the vdW noise. It was found 
experimentally 5 that the vdW force noise in quartz is much higher than in the case of 
metal samples. This difference can be attributed to the dissipation peak around u> t . The 
rough estimate (3.30) of the vdW force noise spectral density in the two slabs model does 
not strongly differ from the experimental data 5 . Some part of this difference may be at- 
tributed to the difference between the simple geometry of the two infinite slabs model and 
the geometry of the real experimental device. 
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APPENDIX: THE MATRIX a mm > 

Let Aj !n i be the coefficient of g n i in the j-th term of the four terms in Eq. (3.11): 1) 
v y + Wy, 2) Vy - w y , 3) v z + w z , 4) v z - w z . Comparing a mml g m ~g* ml and Eq. (3.11) we 
obtain: 



A\ ) ix 

^2, la; 

Al,ly 

A 2 ,ly 

A\,2x 
M,2y 

A s ,i z 

A 3 ,2z 



pqsi 
A 

pqst 



ps\ r 



(s 2 - e 2 p)e ipL + (82 + e 2 p)e- ipL 
(s2~e 2 p)e^ L -(s2 + e 2 p)e-* pL 

(s 2 - e 2 p)e ipL + (s 2 + e 2 p)e- tpL 



ps\ 



(s 2 - e 2 p)e ipL - (s 2 + e 2 p)e 



-ipL 



2pgs 1 s 2 
A ' 

2pS!S% 

A ' 
2 r ps\uu 2 

c 2 A' 6 
2ps 2 uJ 2 

c 2 A' ' 



A2,2x = 
A 2 ,2y = 

~ ipL 7 A 4 

A^2z = 



2p 2 qeis 2 
A 

2e 1 p 2 s 2 2 



2S\S20J 2 



,lz 



-ipL 



c 2 A' 



2s±S2^ 

c 2 A' 



(Al) 



The matrix elements a mm / are presented in terms of A~ m : 



ujuj' I ' c 2 — qq' 



ai x ,i x = A ljlx (q,Lu)A* llx (q' \lu') + 
aiy,i y = A hly (q,u)A* lly (q',u') + 
ai x ,iy = A 1)lx (q,u)Al ly (q',u)') + 
aiy,i x = A hly (q,u)A* llx (q',u') + 
ai x ,2x = A 1Ax (q,u)Al 2x (q',uj') + 
ai y ,2y = A hly (q,Lu)A* h2y (q',Lu') + 
ai x , 2 y = A ljlx (q,uj)A* 12y (q',uj') + 
ai y ,2x = A hly (q,Lu)A* 12x (q',Lu') + 



p(q : uj)p*(q',uj' 

I I 2 / 

lulu J c — qq 



p(q 7 uu)p*(q' 7 uj' 

I I 2 / 

ouou jc — qq 



p(q, Lu)p*(q', lu' 
ujuj' jc 2 — qq' 



p(q 7 u>)p*(q' 7 uj' 

I I 2 / 

lulu jc — qq 



p(q, uu)p*(q' ', uu' 

I I 2 / 

ujuj jc — qq 



p(q 7 uu)p*(q' 7 uj' 
ujuj' /c 2 — qq' 



p(q,uj)p*(q',uj' 
ujuj' /c 2 — qq' 



p(q,uj)p*(q',uj' 
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A 2)lx (q,u)A* 21x (q' ,u'), 
A2,i y (q,Lu)A 2 * Ay (q',Lu'), 
A 2 ,ix(q,uj)A 21y (q',uj'), 
A 2Ay (q,Lu)A^ lx (q / ,LU > ), 
A 2 ,ix(q,uj)A 22x (q',uj'), 
A 2 ,i y (q,Lu)A^ 2y (q',Lu'), 
A 2 ,ix(q,uj)A 22y (q',uj'), 
A 2 ,i y (q,Lu)A 2 ' 2x (q',Lu'). 



(A2) 



au,u = [l - (c 2 /uu')qq']A 3:lz (q, u)A^ lz (q' , u') 
+ (c 2 /uuj')p(q, uj)p*(q', u')A 4ylz (q, u)Al lz (q , u'), 
a 2z ,2z = [l - (c 2 /uu')qq']A 3 ,2z(q,u)Al 2z (q' ,u) 
+ (c 2 /ujuj')p(q, Lu)p*(q', Lo')A 4 ,2 Z (q, u)A^ 2z (q\ u'), 
ai z , 2z = [l - (c 2 /uuj')qq']A 3M (q,uj)A^ 2z (q',uj') 
+ (c 2 /u)U)')p(q, uj)p*(q', u')A^ lz (q, u)A*^ 2z (q', J). 
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